Applied optimal shape design

AbstractThis paper is a short survey of optimal shape design (OSD) for fluids. OSD is an interesting field both mathematically and for industrial applications. Existence, sensitivity, correct discretization are important theoretical issues. Practical implementation issues for airplane designs are critical too.The paper is also a summary of the material covered in our recent book, Applied Optimal Shape Design, Oxford University Press, 2001

Optimum Shape Design Using Automatic Differentiation in Reverse Mode

This paper shows how to use automatic differentiation in reverse mode as a powerful tool in optimization procedures. It is also shown that for aerodynamic applications the gradients have to be as accurate as possible. In particular, the effect of having the exact gradient of he first or second order spatial discretization schemes is presented. We show that the loss of precision in the gradient affects not only the convergence, but also the final shape. Both two and three dimensional configurations of transonic and supersonic flows have been investigated. These cases involve up to several thousand control parameters

An example of a global shock fitting method

AbstractBy means of an example of a single first order equation, we show how a shock can be characterized as the solution of an optimization problem. the optimization problem is solved directly by a gradient calculation and the method of steepest descent. The result is a global shock fitting method which eliminates dispersion. Extensions to the case of a system of equations are discussed

Numerical studies of the Lagrangian approach for reconstruction of the conductivity in a waveguide

We consider an inverse problem of reconstructing the conductivity function in a hyperbolic equation using single space-time domain noisy observations of the solution on the backscattering boundary of the computational domain. We formulate our inverse problem as an optimization problem and use Lagrangian approach to minimize the corresponding Tikhonov functional. We present a theorem of a local strong convexity of our functional and derive error estimates between computed and regularized as well as exact solutions of this functional, correspondingly. In numerical simulations we apply domain decomposition finite element-finite difference method for minimization of the Lagrangian. Our computational study shows efficiency of the proposed method in the reconstruction of the conductivity function in three dimensions

What is the optimal shape of a pipe?

We consider an incompressible fluid in a three-dimensional pipe, following the Navier-Stokes system with classical boundary conditions. We are interested in the following question: is there any optimal shape for the criterion "energy dissipated by the fluid"? Moreover, is the cylinder the optimal shape? We prove that there exists an optimal shape in a reasonable class of admissible domains, but the cylinder is not optimal. For that purpose, we explicit the first order optimality condition, thanks to adjoint state and we prove that it is impossible that the adjoint state be a solution of this over-determined system when the domain is the cylinder. At last, we show some numerical simulations for that problem

Continuum viscoplastic simulation of a granular column collapse on large slopes: μ(I) rheology and lateral wall effects

We simulate here dry granular flows resulting from the collapse of granular columns on an inclined channel (up to 22°) and compare precisely the results with laboratory experiments. Incompressibility is assumed despite the dilatancy observed in the experiments (up to 10%). The 2-D model is based on the so-called μ(I) rheology that induces a Drucker-Prager yield stress and a variable viscosity. A nonlinear Coulomb friction term, representing the friction on the lateral walls of the channel, is added to the model. We demonstrate that this term is crucial to accurately reproduce granular collapses on slopes ≳10°, whereas it remains of little effect on the horizontal slope. Quantitative comparison between the experimental and numerical changes with time of the thickness profiles and front velocity makes it possible to strongly constrain the rheology. In particular, we show that the use of a variable or a constant viscosity does not change significantly the results provided that these viscosities are of the same order. However, only a fine tuning of the constant viscosity (η=1 Pa s) makes it possible to predict the slow propagation phase observed experimentally at large slopes. Finally, we observed that small-scale instabilities develop when refining the mesh (also called ill-posed behavior, characterized in the work of Barker et al. [“Well-posed and ill-posed behaviour of the μ(I)-rheology for granular flow,” J. Fluid Mech. 779, 794–818 (2015)] and in the present work) associated with the mechanical model. The velocity field becomes stratified and the bands of high velocity gradient appear. These model instabilities are not avoided by using variable viscosity models such as the μ(I) rheology. However we show that the velocity range, the static-flowing transition, and the thickness profiles are almost not affected by them

Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows

We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable